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Two-dimensional electron systems with spin-orbit coupling in the proximity of a superconductor 
and a magnetic insulator have recently been considered as promising candidates to realize topological 
superconducting phases. Here we discuss electromagnetic response properties of these systems. 
Breaking of time-reversal symmetry leads to an anomalous Hall effect with a characteristic non- 
monotonic gate voltage dependence and a Hall conductivity that can change a sign as temperature 
is varied. The imaginary part of the Hall conductivity at finite frequency, which shows up for 
example in the Kerr rotation or photoabsorption, can distinguish different topological phases. In 
addition, we demonstrate the existence of magnetoelectric effects associated with the Hall effect; 
in-plane electric fields induce a parallel magnetization and in-plane time-dependent magnetic fields 
induce parallel electric current. 

PACS numbers: 72.25.-b, 85.75.-d,74.78.Fk, 74.78.Na 



I. INTRODUCTION 

A spin-orbit coupling of electrons in solid-state systems 
provides an important link between spin and orbital de- 
grees of freedom, resulting in a wide variety of phenom- 
ena that enable spintronics applications and realizations 
of novel phases of matter. It has been envisioned that 
spin-active components could serve as building blocks of 
future electronic components.^ Electrical manipulation of 
spin, enabled by the spin-orbit coupling, plays a central 
role in these developments. Moreover, spin-orbit effects 
are a crucial ingredient of recently discovered topological 
insulator materials!^ New phenomena and potential for 
applications have attracted enormous interest in topolog- 
ical insulator physics. 

Interesting topological properties exist not only in in- 
sulators but also in so-called topological superconduc- 
tors (TS). Perhaps the most interesting property of TS 
is the existence of Major ana fermions which could serve 
as a platform for topological quantum computation. One 
quantum bit can be encoded in two localized Majorana 
zero modes and computation operations can be carried 
out by braiding such objects. Recently it was proposed 
that topological superconductors can be realized in spin- 
orbit coupled electron systems in the proximity of a su- 
perconductor and a ferromagnetic insulator or in the 
presence of magnetic fields.™^ Besides the interests in 
TS, the interplay of superconductivity, a spin-orbit cou- 
pling and magnetization is interesting in its own right, 
since there already exists a number of experimental real- 
izations coupling a 2DEG with a superconductorP 

In this paper we study unusual electromagnetic 
properties of a Rashba-coupled superconducting two- 
dimensional electron gas (2DEG) in the presence of mag- 
netization perpendicular to the plane. One the most im- 
portant phenomenon resulting from the Rashba coupling 
and time-reversal symmetry breaking due to the mag- 
netization is the anomalous Hall effect. 10 Previously the 



anomalous Hall effect in a spin-orbit coupled supercon- 
ducting system has been studied only numerically in a 
lattice model with a single magnetic impurityP^ Here 
we find that the studied system exhibits a characteristic 
anomalous Hall effect which is non-monotonic as a func- 
tion of chemical potential. The Hall conductivity can 
even change sign when temperature or chemical poten- 
tial are varied. In stark contrast to translationally invari- 
ant chiral p-wave systems which also break time-reversal 
symmetryp^LS the Hall conductivity of the studied sys- 
tem remain finite in the low frequency limit. Moreover, 
the onset frequency of the dissipative part of the ac Hall 
conductivity behaves qualitatively differently in topolog- 
ically trivial and non-trivial phases, enabling an electri- 
cal characterization of phases. Topological nature of TS 
manifests in the quantized thermal conductivity, 16 mea- 
surement of which is challenging in experiments. The 
electromagnetic response studied here provides a signa- 
ture of topological phase transition which should be eas- 
ier to measure. Intriguingly, there exists also magneto- 
electric responses intimately related to the Hall effect; 
in-plane electric fields induce parallel magnetization and 
in-plane time-dependent magnetic fields induce paral- 
lel electric currents. These magnetoelectric effects are 
unique properties of Rashba systems in the presence of 
magnetization and have no counterparts in chiral p-wave 
systems. 



II. MODEL AND ELECTROMAGNETIC 
ACTION 

In this paper we study a 2DEG with a Rashba spin- 
orbit coupling and magnetization in the proximity of a 
s-wave superconductor. The system is described by a 



Bogoliubov-de Gennes Hamiltonian 4 , 

H(k, if) = (e k + a(k x a y - k y a x )) t z + 

M<j z + Acos if r x + Asin ip r y , 



(i) 



where = /i and and Ti are Pauli matrices 

operating in the spin and the particle-hole space, respec- 
tively. Hamiltonian is written in the Nambu basis 

s& (il>kt,i>H,i>-k±i-$-k$) T ' Tne first term corre - 

sponds to the kinetic energy of electrons and holes includ- 
ing the Rashba coupling, the second term is the Zeeman 
splitting due to out-of-plane magnetization and the last 
two terms are proximity-induced superconducting pair- 
ing terms for the order parameter A e lLp . Fundamen- 
tal properties of the model have been discussed in de- 
tail in Refs. [171 and [18] and the effects of disorder have 
been considered in Refs. fl9l and l20l In the absence of su- 
perconducting order, the energy bands are illustrated in 
Fig. [2] for various magnetizations and spin-orbit energies 
Er = a 2 m/2h. We will see below that qualitative fea- 
tures of the Hall response depend heavily on the location 
of chemical potential. The spectrum of ([I]) consists of 
four bands and is symmetric with respect to zero energy 
E = due to the presence of a particle-hole symmetry. 
The positive energy bands 22i(k), £^(k) and the negative 
energy bands i£_i(k), £?_2(k) satisfy E-i(k) — —Ei(k) 
for i = 1,2. The energies are given by 



E 2 1/2 (k) 



el 



2 l.2 



a k 



M 



T2 



NP(4 



A 2 ) 



e 2 k a 2 k 2 



(2) 



with k = |k|, so the spectrum is rotationally symmet- 
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Figure 2: Normal state spectrum as function of the 

in-plane momentum k for different values of 
perpendicular magnetization M. The strength of the 
Rashba coupling is characterized by Er = a 2 m/2h, 
chemical potential is measured from the middle of the 
magnetization induced gap. 



Aq is the scalar potential and A = (A x ,A y ) is the vec- 
tor potential. Below we combine the potentials to a 
single quantity A = (A ,A). It is convenient to per- 
form a gauge transformation through the unitary rota- 
tion U = e l ~^^l^ which leads to the action 

S = - [ d 2 x dr 
2 J 

^ \d r - ieA r z - if (k - eAr z /h, 0) - B • a] *, (3) 



where Aq = Aq 



We have included 



2e ' A ~ A 2e ' 

the Zeeman coupling due to in-plane magnetic fields with 
B = \giJLB(b Xl b y ) 1 where g is the effective ^-factor, \ib 
is the Bohr magneton and (b xi b y ) the in-plane magnetic 
field. Integrating out electronic degrees of freedom yields 
/ V^V^e~ s = e~ s '^ B \ where S'VA, B] is given by the 
determinant of the kernel in Eq. Q. The saddle-point 
expansion up to the second order in the field variables 
produces 

s'[A,b] = l Y J \M-<i)Q { ^)AM+ 



Figure 1: Studied heterostructure formed by 2DEG 
with a Rashba spin-orbit coupling, sandwiched by a 
ferromagnetic insulator and a superconducting layer. 

The coordinates are chosen so that the induced 
magnetization, which is perpendicular to the 2DEG, 
coincides with the z axis and the 2DEG lies in the x — y 

plane. 

fields can be derived by integrating out electronic as well 
as the superconducting phase degrees of freedom, as de- 
tailed in Ref. [12J in the context of a chiral p-wave sys- 
tem. The advantage of this procedure is that it yields 
an action which is explicitly gauge- invariant. The start- 
ing point is the action corresponding to Eq. which is 

S = f cPxdr^ [d T — H(k, cp)] \I//2. The system is cou- calculated for vanishing fields A^ = 0, Bi = and (p = 



B i (-q)Q\ 2 j ) (q)B j + 2B i (-q)Q\ S J(q)A fi (q) 



where we have adopted a convention that repeated in- 
dices should be summed. Quantities Q$ are the stan- 
dard current- current correlation functions describing the 
response to the dressed electromagnetic fields A which 
also contains dynamic of the superconducting phase cp. 
The spin-spin and the spin-current correlation functions 

defined as Q\ j ) (q) = {(Ti<Tj){q) and Q\J{q) = {(JiJ^)(q) 
where the angular brackets stand for two-point functions 



pled to electromagnetic potentials through the substitu- 
tion d T — >• d T — ieAoT z and k — >• k — eAr z /h, where 



and J M denotes the appropriate current operator dis- 
cussed below. The true response is obtained by further 
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integrating out ip, leading to J Vcpe S ^ A ^ B ^ = e s *fd A , B ] 
where the effective electromagnetic action is 

Sefi [A B] = \ [ A^-q)K^{q)A v {q) + 



The off-diagonal response function in the imaginary-time 
representation is 



k.n 



B i (-q)K^(q)B j + 2B i (-q)K%(q)A ll {q) 



(3)/ 
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The kernels are given by 



4^(9) = Q#(?) 



q a qpQ ( ^l{-q)Q [ /3 1 l(q) 
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„(2), v n (2), x 9aqpQiJ{-q)Qpj'(q) 
K\j\q) = Ql/(q) (1) " , 



In the above expression G{iuj n ) = ^ * s the 4 x 
4 Matsubara Green's function of the Hamiltonian = 
i7(k, (/? = 0), (3 is the inverse temperature and ft is the 
area of the system. The summation over i is performed 
over the four energy bands and Pi is a 4 x 4 projection 
operator to the Ei subspace. The trace is evaluated over 
the spin and Nambu indices. Inserting the expressions 
for Green's functions, performing the summation over n 
and analytical continuation to real frequencies yields 



(4) 



Tr [t/ £ Pi J y Pj 
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(rij - rii) , (6) 



where q = (uj,q x ,q y ). The density and current responses 

are given by SJ^q) = ^(^(g) + K$(q)B u (q), . 
where (5J M = (i5p,5J). We have now arrived at purely operators are given by 



where is a Fermi function at energy E{. The projection 

h \ h 2 -eI 



electromagnetic action which is explicitly current con- 
serving q^SJ^q) = 0. 

So far the manipulation has been formal and indepen- 
dent of the detailed form of Eq. Q. Physical properties 
of the electromagnetic response M| are encoded in the 

specific form of functions Q$,. Below we are mainly in- 
terested in the case where the in-plane Zeeman field is 
absent and the response is given by Q$. For super- 
conductors the diagonal responses are non-zero in 
the limit q — >• and uj — >• 0. The spatial components 

at low frequencies are given by the diamagnetic term 

(1) 2 

Q\i ' = , where n s is the supernuid density. The ele- 
ments Qjp give rise to the Meissner effect which screens 
the magnetic field in the bulk. The density-density com- 
ponent Qoo(l) is given by the density-of- states in the 
normal state. In the following we are concentrating on 
effects arising from the off-diagonal spatial components 
Qif 7^ j) which contain signatures from the time- 
reversal symmetry breaking and are responsible for the 
anomalous Hall effect. We are interested in the current 
response in the long-wavelength limit q — » at finite 
uj. According to Q, the off-diagonal physical response 
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functions are given by K, 
K$ = Q^qj/uj (i ^ j) in this limit. 



,(i) K W 



Qifqj/w and 



analogously for P±2 with indices 1 and 2 interchanged. 
With these results the evaluation of the traces is straight- 
forward but tedious. The evaluation is slightly simplified 
by noting that only the terms that involve spin oper- 
ators in Ji contribute to Eq. Also, contributions 
for which i = —j in the sum vanish due to the prop- 
erty (Ei\ J x / y \E_i) = 0. Collecting all the non- vanishing 
terms and converting the /c-summation to integration 
leads to expression 
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(7) 



The expression Eq. ^ is proportional to M and a 2 in- 
dicating that both magnetization and the spin-orbit cou- 
pling are crucial for the off-diagonal term. This also im- 
plies that the Hall effect is independent of the sign of 
spin-orbit coupling a. This is expected since spatial in- 
version takes a — > — a but it does not change the sign 
of the Hall conductivity. The off-diagonal functions sat- 
isfy Q^y (lj) = —Qyx(oj) which is a signature of the Hall 
effect. 



III. OFF-DIAGONAL RESPONSE 



(i) 
y 



In this section we calculate the off-diagonal current- 
current response functions (i ^ j) in the long- 
wavelength limit. The expression for the current oper- 
ator is Ji = + ^ where the antisymmetric 
tensor Sij is defined as en = £22 — 0, £12 = —£21 = 1- 



IV. ANOMALOUS HALL EFFECT 

Here we consider properties of the off-diagonal cur- 
rent response to spatially uniform, slowly varying elec- 
tric fields. As discussed above, the off diagonal re- 
sponse Kij is given by (i ^ j) which leads to the 
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Figure 3: (a): Low- frequency Hall conductivity as a 
function of the chemical potential. The temperature 
dependence of A is assumed to be of the BCS form 

A(T)/Aq = y^l — t^t. (b): Same quantities as in (a) 

but for different parameters. 



Hall effect SJi(uj) = — ^^^ SijEjiuj). For frequencies 
uj <C mmk(E2(k) — Ei(k)), the Hall conductivity is given 
by <r° = lim^o Qxy (w) / 'iuj . Quantity a® y in plotted in 
Figs. |3] (a)- (b) as a function of chemical potential at dif- 
ferent temperatures. In experiments, chemical potential 
could be varied through a gate voltage. For simplicity 
we have assumed that the proximity effect is nearly per- 
fect and the induced superconducting gap obeys a stan- 
dard BCS temperature relation A(T)/Aq = ^1 — 
where A is the gap at zero temperature and T c satisfies 
A = 1.76k B T c . 

The intuition of the behavior of <r!L can be obtained 
by considering the spectrum in the absence of supercon- 
ductivity plotted in Fig. [2] When chemical potential is 
below ~ —En, the Hall conductivity is suppressed since 
the electron density vanishes. As in the normal state 
systems the band curvature effects giving rise to the 
Hall conductivity reach maximum for chemical poten- 
tials located between the two bands in Fig. [2j For larger 
chemical potentials the absolute value of the Hall con- 
ductivity decreases monotonically. The Hall conductiv- 
ity can change significantly even when temperature and 
chemical potential variations are small compared to the 
energy scale of superconductivity A since magnetization 
competes with supeconductivity and suppresses the ex- 
citation gap. Interestingly, the sign of the Hall conduc- 
tivity can change as a function of chemical potential and 
temperature due to the different contributions of band 
curvatures of the four bands. 

The result ^ reduces in the limit A — » and uj — >• 
to the well-known expression of the Berry-curvature 
contribution to the anomalous Hall conductivity of the 
magnetic Rashba model. To see this, first consider the 
first term inside the integral. In the A — >• limit the 
energy bands are given by Ex/2 — \\ e k \ T V M 2 + a 2 k 2 \. 



Using this result one can show that 

e 2 f 1 

o x = —4a 2 M / dkk—— 
h J EiE 2 

( A 2 + e\ E 1 +E 2 \ (ni + n 2 - 1) = 

V Ex+E 2 4 J (E 1 ^E 2 ) 2 

^ M f ^ nA + ^L ( § ) 
2h J Ql (M 2 + a 2 k 2 ) 3 / 2 v y 

where fti denotes part in the k space which satisfy \£k\ < 
V 'M 2 + a 2 k 2 . Using the property n(x)-\-n(—x) = 1 of the 
Fermi function we can write n(£^) +n{E\) — 1 = n(E2) — 
n(—Ei) which further reduces to n(E2) — n{— E\) — 
n(\e k \ + VM 2 + A 2 ) - n(-\e k \ + VM 2 + a 2 k 2 ) in the 
studied parameter regime. Finally, using the property 
n(x\) — n(x2) = n(—X2) — n{—x\), we can drop the abso- 
lute value signs from \e k \ and write n(E2) + n(Ei) — 1 = 
71(62) — n(ei), where £i/ 2 = ^T VM 2 + a 2 k 2 . Thus we 
obtain 

e 2 f a 2 M 
ai = 2h J ai dkk {M i + a^r/i t^ 2 )-n( £l )]. (9) 

Similarly one can show that the second term inside the 
integral of Eq. ^ reduces to ([9| in the complementary 
region in the k space where \£k\ > yfM 2 + a 2 k 2 . Adding 
both contributions together lead to 

e 2 f d 2 k a 2 M r t N t N1 . . 

a " = nJ ^(mwF [n(£2) " n(£l)] ' (10) 

which exactly coincides with the previously known 
expression 10 arising from the Berry curvature in a clean 
system. Thus we see that our general expression for Hall 
conductivity is in perfect agreement with the known ex- 
pression of the normal system. The Hall conductivity is 
a continuous function through the superconducting tran- 
sition. 

The system undergoes a topological phase transition 
between a trivial and a topological phase when /i, A and 
M satisfy the condition /a 2 + A 2 = M 2 P Unfortunately, 
(j® y is not sensitive to the phase transition. The phase 
transition is accompanied by a closing of the energy gap 
in the center of the Brillouin zone such that E'i(O) = 0. 
Vanishing denominator in Eq. ([7]) at the transition point 
is compensated by the vanishing numerator, so a® y is a 
smooth continuous function at the phase transition point. 
Even though the parameters corresponding to Fig. [3] (b) 
are chosen so that the system undergoes a phase transi- 
tion at finite chemical potential, the precise point of the 
transition is not visible. However, we show below that 
the dissipative part of the Hall conductivity at finite fre- 
quency can clearly detect the phase transitions. 

V. IDENTIFYING TOPOLOGICAL PHASES 

At zero temperature the contribution on the third line 
of Eq. ([7]) vanishes, so the imaginary part of the Hall 
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Figure 4: (a): Imaginary part of the Hall conductivity 
as a function of frequency at T = O.lAo, where Ao is 
the gap at T = 0. For all curves M = 1.3A , E R = 2A . 
The threshold frequency for the nontrivial phase (solid 
line) is ujo = 2M above which the imaginary part is 
finite at T = 0. Dashed lines correspond to trivial 
phases. Thermal excitations show up in the positive 
peaks that are more pronounced for trivial phases, (b): 
Zero temperature threshold frequency as a function of 
chemical potential. Existence of a plateau signals a 
nontrivial phase. 



conductivity is finite only for frequencies uj > ujo = 
min/ c (£^i(/c) + E 2 (k)). In the topologically nontrivial 
phase where M 2 — A 2 —/j 2 > 0, the minimum always takes 
place at k = and the threshold frequency is ujo = 2\M\ 
which is independent of the chemical potential. In con- 
trast, in the trivial phase M 2 — A 2 — fi 2 < the threshold 
frequency ujq is always increasing (decreasing) function of 
/x for \i > (/i < 0). In the trivial phase for small chemi- 
cal potentials < A the minimum min^ (Ei (k) + E 2 (k)) 
also takes place at k = but the threshold frequency is 
ujo = 2y // A 2 + /J?. Thus the nontrivial phase is charac- 
terized by horizontal plateaus in the plot ujq vs. /i, termi- 
nating at critical values /i c = ±VM 2 — A 2 corresponding 
to phase transitions to the trivial phase, as illustrated in 
Fig. [4j By measuring the onset frequency ujo while vary- 
ing chemical potential, it is possible determine whether 
the system is in the topological or in the trivial phase. 

Experimentally it has been verified that the proxim- 
ity induced gap in 2DEG may be at least of the order 
of Ao ~ 0.1 meVP Assuming that M can be made also 
be of that order (and somewhat larger to enter the topo- 
logical phase), the threshold frequency ujo is of the order 
of few tens of GHz. This provides the upper limit of 
the frequencies of the interest in the identification of the 
topological phases. 

Characterization of TS phases by electrical means like 
proposed here or in tunneling experiments proposed in 
Ref PH is important since measurement of quantized ther- 
mal conductivity, the natural topological invariant of TS, 
is very difficult. Both the real and imaginary part of 
the Hall conductivity can be probed optically through 
the Kerr effect, which has been employed to characterize 
broken ti me-re versal symmetry in the p-wave candidate 
SrsRuOztP^This technique requires comparing polar- 
izations of incident and reflected electromagnetic waves. 



VI. MAGNETOELECTRIC EFFECTS 

An interesting consequence of the non-zero Hall con- 
ductivity is the existence of unusual magnetoelectric ef- 
fects. First we consider the response of the system 
to in-plane electric fields Ei(uj) = iujAi(uj), which re- 
sults in magnetization parallel to the applied fields. To 
derive this effect, we consider induced magnetization 
pSi = f"Q which can be expressed in the standard 
linear-response theory as 



(ii) 



where K^\uj) is the spin-current response function (jij). 
In the long wavelength limit q — » at finite uj the re- 



sponse function reduces to K^\uj) = Q\j } (uj) with the 
imaginary-time representation given by 

Q[f(iuj m ) = -^—^2TY[o-iG(i(uj m + v n ))JjG(iv n )). 



2np 



(12) 



Since the off-diagonal current response ([5| arises solely 
from the spin part of the current operator, we see that 
the diagonal component of expression ( [l2| ) is give by 

Qyyi 00 ) ~ • Therefore the magnetization parallel 

to the electric field is related to the Hall conductivity as 



h 2 iQiWuj) h 2 o 

= 2^^^^ H = 2^^ H ' (13) 



PSy M 



where we took the low-frequency limit. Since the sys- 
tem is rotationally invariant, similar relation holds also 
in the x direction. This effect can be intuitively under- 
stood as follows. The application of electric field induces 
Hall current in the perpendicular direction. Because of 
Rashba coupling, current is accompanied by magnetiza- 
tion perpendicular to current! 2 -^ Thus, the application of 
an electric field results in magnetization parallel to it. 
This phenomenon is a consequence of the Rashba cou- 
pling and a finite out-of-plane magnetization M ^ 0. 
No analogous effect exists in chiral p-wave superconduc- 
tors. The Rashba coupling also results in a previously 
discovered magnetization perpendicular to applied elec- 
tric fields which remains finite also for vanishing out-of- 
plane magnetization M = 0.^ 

There exists also a Zeeman-type magnetoelectric effect 
closely related to the Hall effect. Suppose that the system 
is exposed to an in-plane magnetic field B x which couples 
to the spin of the particles through a Zeeman term B x o x . 
Similarly as before, we can analyse the linear response of 
the current to the Zeeman field. Analogously we find 

SMw) = K^MBjiu,), (14) 

Therefore, for similar reasons than discussed above, the 
relevant response function can be expressed in the long 
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wavelength limit as K^ 3 J (uj) = Q^J(uj) = 



Thus the parallel part of the current response to the ap- 
plied Zeeman field is given by 



5J x (uj) 



iwB x {uj) 



h 



a xy d t B x (uj), 



(15) 



where the last form is valid for low frequencies. The mag- 
nitude of the effect depends on the Hall conductivity and 
the ratio of the magnetic moment and the spin-orbit con- 
stant a. Interestingly, this magnetoelectric response flips 
the sign upon the changing of the sign of a, in contrast to 
the Hall effect which is insensitive to the sign. It should 
be noted that the parallel current (15) can be distin- 



guished from the field induced Meissner current and the 
previously studied paramagnetic current!^ by its different 
directional and functional dependence on the magnetic 
field. For a linearly polarized in-plane magnetic field cur- 
rent (15) is parallel (or antiparallel) to the applied field, 



not perpendicular to it like the other contributions. In 



further contrast, for static magnetic fields current (15) 
vanish. 



VII. CONCLUSIONS 



We studied electromagnetic response properties of su- 
perconducting two-dimensional electron systems with 
Rashba spin-orbit coupling and perpendicular magneti- 
zation. We focused on the anomalous Hall effect and 
closely related magnetoelectric effects. The anomalous 
Hall effect is finite in the low-frequency limit, exhibit- 
ing a non-monotonic behavior as a function of chemical 
potential. Moreover, the frequency-dependent Hall con- 
ductivity enables a purely electrical characterization of 
different topological phases. We found two magnetoelec- 
tric effects directly related to the Hall effect, one leading 
to a parallel magnetization as a response to in-plane elec- 
tric fields and the other giving rise to electric currents 
parallel to time-dependent in-plane magnetic fields. 
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